Plasticity and reversibility of structural transitions in a model solid 



Arya Paul, Surajit Sengupta 
S.N. Bose National Centre for Basic Sciences, Block JD, Sector III, Salt Lake, Calcutta 700 098, India 

Madan Rao 

Raman Research Institute, C. V. Raman Avenue, Bangalore 560 080, India 
National Centre for Biological Sciences (TIFR), Bellary Road, Bangalore 560 065, India 

We formulate a phenomenological elasto-plastic theory to describe a solid undergoing a structural 
transition from a square (p4mm) to an oblique (p2) lattice in two dimensions. Within our the- 
ory, the components of the strain may be decomposed additively into separate elastic and plastic 
contributions. The plastic strain, produced when the local stress crosses a threshold, is governed 
by a phenomenological equation of motion. We investigate the dynamics of shape of an initially 
square solid as it is cycled through a transformation protocol consisting of (1) a quench across the 
transition (2) deformation by an external stress and finally (3) reverse transformation back to the 
parent state. We show that shape recovery at the end of this cycle depends on crucially on the 
presence of plasticity in components of the strain responsible for the transformation. 

PACS numbers: 



I. INTRODUCTION 

The question of the degree and efficiency of shape 
reversibility of solids undergoing reversible Martensitic 
transitions [1] is clearly of considerable technological im- 
portance considering the diverse applications of shape 
memory or "smart" materials [2]. 

The necessary criterion for a material in order to ex- 
hibit shape-memory is that the symmetry group of the 
product phase forms a sub-group of the parent and the 
transformation leads to symmetry breaking [3]. In a re- 
cent publication we had explored the possible sufficient 
conditions for shape reversibility using a model system 
in two dimensions as a test case[4(. In our model sys- 
tem it was possible to influence the nature of the trans- 
formation by tuning the parameters of the interaction 
potential We found that, while a square (p4mm) to a 
rhombic (in general, oblique - p2) lattice was reversible [6], 
the transformation to the more symmetric, triangular 
(p6) lattice was not, in accordance with the conclusions 
of Ref.[3]. However, we found that the question of re- 
versibility was linked to the extent and nature of plastic 
or non-affine zones (NAZs) present within the transform- 
ing crystal. The positions of particles which belong to 
the NAZs cannot be obtained from a simple afflne (shear 
and/or scaling) transformation of the parent lattice. We 
identified two kinds of NAZs associated either with the 
non-order parameter (NOP) (e.g. volumetric) strain or 
with the order-paramter (OP) shear strain which trans- 
forms the square to a rhombic lattice. We argued that 
the p4mm ^ p2 produces only NAZs of the first kind 
while during a p6 — >• pAmm transition, we gave indirect 
evidence that both kinds of NAZs are likely to be pro- 
duced leading to irreversibility. We concluded therefore 
that for reversibility it is sufficient that NAZs associated 
with OP strains not be produced. To support our claim, 
we showed that eliminating NAZs by making the solid 
stronger or the system size smaller than the typical dis- 



tance between NAZs makes even the p4 — >• p6 transfor- 
mation reversible in apparent contradiction to Ref.[3]. 

In this paper, we lend further support to this claim 
by investigating this problem from the point of view of 
a phenomenological elasto-plastic theory [7] of the square 
to oblique transformation. The elasto-plastic theory has 
been shown earlier to reproduce the qualitative features 
of the results of our MD simulations [5]. One of the ad- 
vantages of this approach is that plasticity in the NOP 
or the OP sector may be introduced systematically and 
their effects observed explicitly. 

Our main results are as follows. We study the re- 
versibility of shape of a solid which is first quenched from 
a square into the rhombic phase, deformed by an external 
stress and finally brought back to the square phase (see 
Fig.l). The external control parameters are the temper- 
ature and stress. We compare the initial and final shapes 
of the solid and look for congruence. We show that the 
presence of plasticity exclusively in the NOP sector does 
not affect shape reversibility while the presence of plas- 
ticity in the OP sector triggers irreversibility in shape 
changes. 

The next section is dedicated to a quick review of our 
elasto-plastic theory where we discuss the formalism and 
how to deduce shape at each instant of time from the 
knowledge of strains at that instant. In section III we 
present our results and we conclude in section IV indi- 
cating possible future directions. 



II. ELASTO-PLASTIC THEORY 

In this section we formulate our elasto-plastic con- 
tinuum theory for structural transitions in solids. Our 
theory allows us to determine the microstructure and 
the overall shape of a solid undergoing a square to a 
general oblique lattice. Unlike previous work in this 
subject [!| we allow for the presence of non-affine, or plas- 
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tic deformations [El 0- 

Our main assumption in what follows is that the com- 
ponents of the strain tensor [8| may be additively de- 
composed into elastic (affine) and plastic parts [10], viz. 
e ij = e tj + e fj wnere hj — X ,V- The plastic strains 
efj represents the total contribution from, in general, 
space and time dependent defect fields, which, in the 
spirit of our continuum approach, need not be resolved 
into individual, microscopic defects. Nevertheless, de- 
fects introduce multivaluedness in the particle displace- 
ment field u = (u x (r),u y (r)) and cause the elastic parts 
of the strain to violate the St. Venant's compatibility 
condition Mil, 
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where e A and e p represents the affine (elastic) and the 
plastic strain tensors respectively. Note that Eq.(pQ) im- 
plies that the total strain e^- does satisfy the compatibil- 
ity condition. 

The structural transition is driven by the non-linear 
response of the solid to one or more components of the 
strain- the order parameter (OP). For example, for the 
(p4mm — >• p2) transition, the order parameter strains are 
e 2 = e xx — e yy and e 3 = e xy = e yx while the remaining 
volumetric strain e\ = e xx + e yy is a non-order parameter 
(NOP) strain[6|. We shall first consider the case when 
plasticity exists only in the NOP sector. 

The transition is described by the following free energy 
functional: 



T = X - J dxdy [ai(ef) 2 + a 2 e 2 + a 3 e 2 3 - b 3 e\ + d 3 $) 

+ Cl (Ve^) 2 + c 2 (Ve 2 ) 2 + c 3 (Ve 3 ) 2 ] . 

Note that in Eq.(|3]) we have made the simplifying as- 
sumption that the product oblique lattice is actually a 
rhombus so that the equilibrium value of e 2 =0. This 
choice is motivated by our MD simulation of a particular 
model solid [5] though we do not doubt that our theory 
may be easily extended to the general case. The three 
elastic constants ai,a 2 and a 3 define the linear elastic- 
ity of the square (parent) phase. The connection with 
external control parameters such as temperature (T) is, 
as usual, through the temperature dependence of these 
coefficients especially a 3 ex (T - T*)/T* where T* is the 
limit of stability of the square solid. Reducing a 3 by cool- 
ing the solid stabilizes the rhombic phase (see Fig. 1). 
The rest of the coefficients parametrize non-linearities 
and may be taken to be constants. 
The Lagrangian is given by [9], 
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To obtain the equation of motion in the displacement 
fields, we need to solve the Euler Lagrange equation: 
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The Rayleigh dissipation functional [H 1Z is given by 1Z — 
- j [£(e^) 2 + fte 2 , + 7e 2 ] dxdy, where the coefficients £, k 
and 7 are the corresponding viscosity coefficients of the 
system. 

Using Eqs.®, (j3]) and (J4j), we obtain the following 
equations of motion for the affine strains ef , e 2 and e 3 : 
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operator W 2 is defined as W 
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d 2 /dx 2 — d 2 /dy 2 . Since there is symmetry breaking in 
the space of the OP strains, the dynamics of e 2 and e 3 
- the broken symmetry variables - is much slower than 
that of the affine NOP strain ef. Therefore ef reaches a 
steady state much faster compared to e 2 and e 3 and we 
can simplify Eqs.(|5]l6]) by assuming that is slaved, at 
all times, to the OP strains e 2 and e 3 . Also note that, 
the St. Venant's compatibility condition, which in two 
dimensions can be explicitly written as 



V 2 ei - W 2 e 2 - 4V x V v e 3 = 0, 
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may be used to eliminate e 2 from Eqs.([5][6]). 

The final set of equations which therefore has to be 
solved to obtain the affine strains are 
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Lastly, the dynamics of the plastic strain ef needs to 
be specified. We assume that this is given by the follow- 
ing simple phenomenological equation of motion used by 
us (Fig. 1 inset) is, 
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where the local volumetric stress o\ = ST/Sef, and we 
have chosen a Newtonian ansatz with a simple threshold 
criterion with yield stress o~\ c for simplicity. Note that 
EqfTTl also includes strain diffusion. 

The continuum elasto-plastic theory described above 
has been successfully used in the past[7[ to describe 
microstructure selection following a quench from the 
square to the oblique phase. For large values of h\ in 
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the affine strains, obtained from this Lagrangian are 



FIG. 1: The homogeneous part of the free energy func- 
tional T{ez) as a function of the OP strain for differ- 
ent values of the parameters as and the external stress (73 
viz. CL3 = 0.5,(73 = (I), CL3 — 0.1,(73 = (II) and 
a 3 = 0.1,(73 = .1 (III)- The rest of the parameters 63 = — 1 
and d% = 1 throughout. Inset shows a schematic representa- 
tion of the typical dynamics of plastic strain ef as a function 
of the local stress a used in our calculations. Similar dynam- 
ics is assumed for volumetric (NOP, i — 1) as well as shear 
(OP, i = 3) plasticity. The threshold stress <Ji C is, of course, 
different in the two cases. 



Eq. ([TT]) ef is produced only at the parent-product inter- 
face within localized and transient NAZs. Once coarse- 
grained over distances large compared to the typical size 
of these NAZs, the twinned microstructure obtained in 
this limit become identical to that seen in "strain-only" 
approaches [|j which do not involve plasticity. On the 
other hand, for small values of hi, ef proliferates the 
entire solid and effectively screen non-local elastic in- 
teractions, completely changing the microstructure from 
twinned martensite to un- twinned ferrite. In this limit 
the coarse- graining length grows to become of the order 
of the system size and such microstructures cannot be 
described within the earlier strain only theory 0|. These 
results closely mimic those observed in our MD simula- 
tions of the square ±= oblique transformation in a model 
solid @. 

In this paper we wish to study shape recovery follow- 
ing a quench- deform- reheat transformation cycle for a 
solid undergoing a structural transition in the presence 
of plasticity. We therefore need to include the effect of 
external stress by using a modified Lagrangian, namely, 



F+fxUx + fyUy (12) 



where f x and f y are the components of the tangential 
forces along the x and y direction respectively applied at 
the edges of the system. The dynamical equations for 
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Note that, the external shear stress <r 3 = + ^7 does 
not affect the steady state equation for the affine NOP 
strain ef. The phenomenological equation for the plastic 
strain also remains same as well. 

If, in addition, plasticity is associated with the OP 
strain, cr 3 may also give rise to plastic flow in the solid. 
In this work, we consider this possibility as well. In order 
to incorporate plasticity in the OP sector, we proceed in 



an identical manner. We write e 3 



- eg . We replace 



e 3 with ef in Eq.(|9]) and Eq.([T3]). the dynamics of e\ 
is assumed to be given by a phenomenological equation 
similar to eqn ([TT]) . Finally, as before, the total strains e\ 
and e 3 satisfy the St. Venant's compatibility condition 
Eq.®. 

Having specified the dynamics of the strains, it is nec- 
essary to be able to compute the shape of the solid from 
the strain fields in order to investigate shape reversibil- 
ity. For this, we make use of the Kirchoff-Cessaro-Voltera 
|ll| relation 

u(r) = / [E(l) + (1 - r) x Vi x E(l)] • d\ (14) 

where E(l) is the strain tensor at position 1. The line inte- 
gral is along any arbitrary contour C(i*o,r) from a fixed 
point of the deformation 1*0 to the point of interest r. 
The displacements, so evaluated, are valid upto a global 
translation and global rotation, that can be viewed as in- 
tegration constants. In two-dimensions, Eq. ([T4j) reduces 
to, 
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Evaluating the displacements u x and u y at the bound- 
aries will generate the shape of the system subjected to a 
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FIG. 2: (a)(i) Plot of showing the quenched microstruc- 
ture. The parameters are p — 1,<23 — 0.04, 63 — — 6, cfe — 
4, c 3 = 4, 01 = 4, ci 0.04, a 2 0.04, c p 0.04, 7 = 20, £ = 
0.4, aic = 1 and hi = 10. (ii) Corresponding plot of the affine 
strain ef. (b)(i) The micro- structure and (ii) ef when the 
system is placed under an external shear stress. The param- 
eter (73 = 0.5. (c)(i) The order parameter when the system is 
"heated" by increasing a 3 from 0.04 to 4. The entire system 
goes into a square phase phase with ez = 0. (ii) The remnant 
volumetric strain ef averages to zero. The color range is from 
-1 (black) to (brown) to 1 (yellow) for (i)(a),(b)&(c). For 
(ii)(a),(b)&(c) the corresponding range is from -.1 (black) to 
.1 (yellow). 

global rotation and a global translation. Note that total 
strains, which satisfy the compatibility constraint Eq.flH), 
are used to evaluate the displacements. 

III. RESULTS 

As mentioned before, we look for shape reversibility at 
the conclusion of the following transformation protocol. 
Firstly, a finite size solid initially in the shape of a square 
with a parent (square or pAmm) crystal structure (Fig. 1 
curve I) is quenched to get a rhombic martensitic product 
with a p2 structure (Fig. 1 curve II). Next, the product 
phase is deformed (Fig. 1 curve III) with the help of an 
external shear strain 0-3. Finally, the deformed system is 
transformed back to the parent square crystal structure 
(Fig. 1 curve I). All throughout we keep track of the 



shape using Eqs. (fT5|) and (fTBj) . If the system returns 
to the square shape after heating, we conclude that the 
shape is reversible, otherwise not. We have performed 
this transformation protocol for two cases: Case I- when 
there is only NOP plasticity and no OP plasticity Case 
II- when plasticity is present in both NOP and OP sector. 
We give below these results for the two cases, one after 
the other. 



A. Case-I 

We first consider the case where plasticity exists only 
in the NOP sector. We need to solve Eqs.® and ([13)) to 
get the desired quenched microstructure from the parent 
square phase. Our initial values of es(r) and ef(r) are 
zero everywhere. For the plastic strain ef , we assume 
ef (r,0) = S(x-L/2)(y-L/2) at t = 0. When quenched 
below the transformation temperature, the OP strain es 
shows the formation of twins, originating from the center 
of the simulation box where the initial value of the plastic 
strain was nonzero [see Fig(2)(a)(i)]. Fig(2)(a)(ii) shows 
the affine NOP strain ef. The calculated shape of the 
quenched solid is shown in Fig(3)(a)(i). Although alter- 
nating twins are present, the overall shape of the solid 
continues, on an average, to be a square. 

We then load the microstructure using an external 
shear stress so that only one of the product variants is 
favored. The system prefers to go to a single variant in 
order to lower the free energy. We stop the evolution af- 
ter the system has undergone sufficient deformation. The 
plots of the OP strain es and the affine NOP strain ef is 
shown in Fig(2)(b). The shape of this deformed system 
is shown in Fig(3)(b)(i). 

Having deformed the system, we transform it back to 
the parent phase by increasing as (temprature). The 
plots of the OP strain es and the affine NOP strain ef 
is shown in Fig(2)(c). The whole system relaxes to zero 
OP strain. The final shape of the system is plotted in 
Fig(3)(c)(i). The system goes back to the parent phase, 
i.e., a square phase, and, at the same time recovers its 
shape at the end of the complete transformation cycle. 
The plasticity in the NOP sector does not affect shape 
recovery because the ef alternate in sign and average to 
zero over the entire sample. The shape is determined 
mainly by the OP strain which is completely reversible. 



B. Case II 

We now look into a system undergoing square to rhom- 
bic transition including plasticity in the OP sector as 
well. As before, we start with a solid in the shape of a 
square and in the p4mm phase with a delta function NOP 
plasticity at the center of the simulation box and let the 
system evolve without any external stress to obtain the 
microstructures. The plots for the affine OP strain ef 
and the affine NOP ef are shown in Fig. 4(a). Note that 
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FIG. 3: (a)(i) - (c)(i) The shape of the system corresponding 
to quenched, loaded and heated microstructure in Fig. (2) 
and (a)(ii) - (c)(ii) Fig. (4). 



these, as well as the calculated external shape, shown in 
Fig. 3. (a)(ii), are not significantly different from those 
obtained in Case I. 

The quenched microstructure, so obtained, is then, as 
before, deformed with an external shear stress as in the 
previous case. Unlike Case I, however, the external de- 
formation now has two effects. Firstly, as before, the ex- 
ternal stress favors one of the twin variants over the other 
and therefore alters the microstructure by changing the 
distribution of the variants. Secondly, the external stress 
may cause the solid to flow, so that a part of the de- 
formation, would now be plastic. The amount of strain 
which appears as the plastic strain in the OP sector, of 
course, depends on the parameters eg as c - the threshold 
stress and hs - the rate of plastic strain production. In 
the present work, we have used a low value of the stress 
threshold to over-emphasize the plastic contribution for 
the sake of illustration. In reality, the contribution of the 
plastic strain may be smaller. The plots of the affine OP 
strain ef and the affine NOP ef are shown in Fig(4)(b). 
We use the obtained strains to evaluate the displacements 
at the boundary and hence the shape of the system which 
is plotted in Fig(3)(b)(ii). 

The system is then transformed by increasing as and 
the resulting ef and ef are plotted in Fig(4)(c). There 
are remnant (plastic) shear strains in the system even 
after the solid is fully transformed so that the system 
fails to go back to the initial square shape as can be 
seen from the plot of Fig(3)(c)(iii) although the crystal 
structure is p4mm everywhere as expected. The shape 
of the solid is now not only determined by the affine OP 
strain ef but also, by the plastic part ef which is not 
recovered on increasing as. Ultimately this contribution 
causes shape irreversibility. 
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FIG. 4: (a)(i) Plot of ef for the quenched microstructure. (ii) 
Corresponding plot for ef . (b)(i) Plot for ef for the deformed 
microstructure. (ii) Corresponding plot for ef . (b)(i) Plot for 
ef after heating the system, (ii) Corresponding plot for ef . 
(T3c = 2 and /13 = 1. All other parameters are same as that 
of Fig(2). For (i)(a)-(c), the color range in from -1.5 (black) 
to 1.5 (yellow) For (ii)(a)-(c), the color range goes from - 
0.1 (black) to 0.1 (yellow). 



IV. SUMMARY AND CONCLUSION 

Using atomistic computer simulations of a model solid, 
we had shown in Ref.[4] that group-nonsubgroup trans- 
formations are irreversible because plasticity is automati- 
cally generated in the OP sector during such transforma- 
tions in agreement with the conclusions of Bhattacharya 
et a/.[3|]. We showed further that if OP parameter plas- 
ticity is suppressed by increasing the yeild stress or in a 
small system, even group-nonsubgroup transformations 
may become reversible. In this paper, we complete this 
line of reasoning by showing explicitly that for shape re- 
covery it is sufficient that all plastic strains be associated 
exclusively with the NOP sector of the transformation. 
If, on the other hand, OP plasticity is present, complete 
shape recovery is impossible. We show this here using an 
elastoplastic theory where it is possible to control the de- 
tailed nature and extent of the plastic strains by varying 
parameters. If OP plasticity is therefore "artificially" 
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introduced in a group - subgroup (p4mm ^ p2) the 
normally reversible transformation becomes irreversible. 
Our results may be verified experimentally using shape 
memory alloys with known and non-trivial yeild criteria 
such that plasticity in either sector may be controlled 
independently of each other. We await such systematic 
studies on the effect of plasticity on shape reversibility. 

Before we end, we would like to discuss several possible 
extensions of this work which we take up one by one as 
follows. 

Stress induced martensite: Within our formalism it is 
straightforward to consider the effect of plastic defor- 
mation on stress induced martensitic transformations [lj. 
Stress- strain hysteresis, which can be routinely measured 
in such systems may show nontrivial effects due to plas- 
ticity at various spatial and temporal scales. The sepa- 
rate effects of NOP and OP plasticity should be explicit 
in such studies which are planned in the near future. 

Scale dependent reversibility: Reversibility of a trans- 
formation is, obviously, scale dependent. We have stud- 
ied here is the question of shape reversibility, namely, 
reversibility at the largest scale available to the system. 
In contrast, complete microscopic reversibility would im- 
ply that the very positions of atoms be recovered as the 
transformation is reversed - an impossibility considering 



the identity of atoms and thermal noise. At intermediate 
scales, one may ask whether microstructural features are 
reversible or not below a certain "irreversibility length" . 
Such calculations are in progress and will be published 
elsewhere. 

Complex dynamics and defect reorganization: The non 
uniform plastic strain fields created during the transfor- 
mation are redistributed as the solid is cycled through the 
transformation protocol many times. It is legitimate to 
ask whether this distribution has a steady state and what 
its dynamic properties may be. Defect redistribution and 
aging during transformation with nontrivial d yna mical 
signatures have been observed in real systems[12]. We 
plan to study such questions, as well, within our elasto- 
plastic approach in the future. 
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